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Recent CLEO results on the observation of the T(1D) state and a search for the ^(lS) state 
from a data sample of 4.7 x 10 6 T(3S) decays collected with the CLEO-III detector, and a 
measurement of the rj production spectrum using 1.9 x 10 6 Y(IS') decays collected with the 
CLEO-II detector are described. 



1 Introduction 



CLEO is concluding a 1-year program to collect ~ 4fb~ 1 of integrated luminosity distributed 
over the T(15),T(2S) and T(35) resonances. These datasets will augment the world supply of 
bb quark-antiquark bound state data for these resonances by more than an order of magnitude. 
The data is produced by the symmetric e + e" collider CESR and collected by the CLEO-III 
detector, a configuration of the CLEO detector with excellent electromagnetic calorimetry and 
enhanced particle identification capabilities. 



2 T(1L>) Observation 

Heavy quarkonia (e.g., bb quark pairs) are an appealing physical system for probing the strong 
interactions. Bound states below open flavor threshold are narrow, leading to long lifetimes and 
little mixing between excitation levels. The constituent quarks are essentially non-relativistic, 
making the bound state system particularly appropriate for analysis by lattice QCD calculations 
and by effective theories via potential models. 

We present evidence for the observation of the T(1D) state, the first observation of a narrow 
bb bound state since the early 1980's, using the following radiative 5-stage cascade decay that 
ultimately produces 4 photons and a lepton pair in the final state: 



T(35) - 7Xfe(2Pj) 

Xfe (2Pj) - 7 T(LD) 

T(1L>) -> TXo(l^) 

X6 (lPj) - 7 T(LS) 

T(15) -» e+e",^" (1) 

The theoretical product branching fraction 1 for this sequential decay is 3.76 x 10~ 5 . The 
dominant background reaction is T(3S I ) — > 7r°7r°T(lS'); T(15) — > e + e~,fi + fi~. An additional 
background mode, the same cascade chain as Eq. 1 except that the intermediate T(1D) is 
replaced by the T(25) state, must also be rejected. 

Measuring the branching fraction for the background T(35) — > 7r°7r°Y(lS) transitions is a 
useful check of the technique used to reconstruct our cascade signal mode. A key analysis variable 
X^o^o used to reconstruct T(3S) — > 7 77 7Y(1S') decays is the minimal 7r° mass deviation chi- 
squared formed from among the three possible photon pairs formed from four individual photons. 
After cutting on X^o^o, the ratio AM/a(AM) is formed from the deviation of the four-photon 
recoil mass from the T(1S) mass AM = M recoi u^ — Myr 1S \ and its expected mass resolution 
cx(AM). The distribution AM/a{AM) is shown for data in Fig. 1 fit with a Monte Carlo 
determined shape. Combining the yield of 737±28 events with a 13.6% Monte Carlo determined 
efficiency implies a branching fraction B(T(35) -> 7r°7r°T(lS)) = (2.33±0.09±0.16)%, consistent 
with, but more precise than, previous measurements 2,3 . 
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Figure 1: Distribution of the recoil mass deviation from the T(3S) mass for T(3S) — » 7r°7r°T(lS) candidate 
events. Dots are data, solid line is fit. 

Confident that our technique for reconstructing multiple photons in a single event is robust, 
we select T(1D) events by constructing a chi-squared based variable x 2 (^D) that is a function 
of the unknown T(1D) mass and the total angular momentum J of the \b states on either 
side of the T(1D) in the cascade chain. Other chi-squared based variables are used to reject 
background. Surviving events are shown in the left hand side of Fig 2 as the histogram with the 
signal region x 2 (l-D) < 10. The fit of the signal on top of the background is shown as the heavy 
and dashed lines, respectively, producing a significance of 9.7a. 

The Y(1D) mass is estimated by plotting for each event in the signal region of the left hand 
plot in Fig. 2 the value of the mass that minimizes x 2 (lD). Monte Carlo studies predict a 
peak at the true mass plus two satellite peaks due to mis-reconstructed photons from shower 




Figure 2: The left hand side is a fit to the Xid f° r data. The solid line is a fit to the signal contribution on top of 
backgrounds. The right hand side is a fit of a two-peak structure to the most-likely mass taken from the signal 
region of the left hand plot. 

leakage in the calorimeter. The mass distribution for data is shown as the histogram in the 
right hand plot of Fig. 2. Although the T(1D) can be produced in any of the J = 1, 2, 3 states, 
our selection efficiency is largest for the T(1D±) and T(1Z>2) states while theory predicts that 
the T(1Z?2) state should be preferentially produced in the decay cascade. Fitting the histogram 
with two peaks, each with two satellites, produces a peak at Mi ow = 10161.2 ± 0.7MeV/c 2 and 
at M high = 10174.2 ± 1.3MeV/c 2 , with an overall 58% C.L. Fitting under the assumption that 
there is no low mass peak produces a 0.04% C.L. and the difference in significance between these 
two fits implies a significance of the low mass peak of 6.8<r. Fitting the histogram with just a 
single peak plus two satellites yields a 43% C.L. fit with a peak at 10162.0 ± 0.5MeV/c 2 . 

Averaging over the different fits yields a T(LD) mass of 10162.2 ± 1.6MeV/c 2 with the 
J = 2 state favored, but the J = 1 state not excluded. Based on our measured T(1D) yields 
and Monte Carlo determined efficiency, the preliminary overall branching fraction for the 5-stage 
decay cascade of the T(3S) via the T(LD) is (3.3 ± 0.6 ± 0.5) x 10" 5 . 

3 Search for r] b (lS) 

The spin-singlet state of the bb system, including the bb ground state r]b(lS), has yet to be 
observed. The reaction e + e~ — > 7* — > bb produces bb states with the same quantum 
numbers as the photon so the ^(l 1 ^) state cannot be produced directly at CESR. However, 
the state is accessible through decays from the spin-triplet T(nS) states via magnetic dipolc 
(Ml) radiative transitions which occur between states with opposite quark spin configuration 
but the same orbital angular momentum. The rate for Ml transitions between the T and the 
r)t, is given by 4 

r[T(nS) ->7ffo(n'S)] oc I 2 ^, (2) 

where n and n' are the principal quantum numbers of the respective bb systems, I is an overlap 
integral between initial and final bb states and is the energy of the emitted photon. 

With a data sample of 4.7 x 10 6 Y(3<S) events, CLEO searches for transitions with n 7^ n 1 to 
avoid photon energies in the 100 MeV range which are susceptible to background contamination 
from 7T° decays. The inclusive photon search technique for T(35) — > 7% (IS) is tuned using the 



electric dipole (El) transitions Xb(2P) — > T(15) where the photon energy is ~ 800 MeV, similar 
to what is expected 4 for the signal Ml transition. 

The inclusive photon spectrum between approximately 550 MeV and 1 GeV is shown in 
the left hand side of Fig. 3 where the dots indicate data and the dashed line is background. A 
logarithmic energy scale is used to preserve energy resolution. The central peak is a superposition 
of three gaussians for the three Xft(2-Pj) El transitions and the search window is also indicated. 




Figure 3: The inclusive photon spectrum in the El peak and the r/i, search region before (left) and after (right) 
background subtraction. 

The right hand side of Fig. 3 shows the background subtracted central peak and the search 
region fitted with a Crystal Ball lineshape. Multiple fits are performed for various photon 
energies in the search region and no signal is seen. This lack of signal is converted into a 90% 
C.L. upper limit shown in Fig. 4, excluding a variety of phenomenological models 4 . 



Hyperfine Splitting: M Y(1S) - M^,^ (MeV/c ) 
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Figure 4: Upper limits on S(T(3S) — ► ^yrjb(lS)) as a function of E 7 . Model predictions are from Ref. [4]. 



4 T(15) -> r/X Production 



An unexpectedly large 1? — > r/'X^ inclusive decay rate in the momentum range 2 < < 2.7 
GeV/c observed by CLEO 5 and BABAR 6 is possibly described by an anomalously large coupling 
between the rj and two gluons 7 ' 8,9 in the underlying b — » sg mechanism thought responsible for 
rj production. Here the process b — > sg 1 is followed by the two gluon coupling to the rj shown in 
Fig. 5. 
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Figure 5: Feynman diagram for b — > s<?r/. 

The effective rf g*g coupling can be written as 10 

Htfy^fkP^Z, (3) 

where H(q 2 ) is the r)'g*g transition form factor, q = Pb~ Ps is the four-momentum of the virtual 
hard gluon (g*) and k is the momentum of the soft gluon (g). Workers 10 have shown that the 
q 2 region relevant to the process b — > sgrj' is accessible in high energy rj production in T(IS') 
decays so that constraints can be placed on H(q 2 ) from the fast rj spectrum in T(1S) — > ggg 
decays. Various choices 8,10 are available for H(q 2 ). 

We use 1.86 x 10 6 T(LS') decays collected with the old CLEO II detector and detect rj mesons 
through the decay channel rj — > r/7r + 7r~(r/ — ► 77). This data sample can be divided into three 
parts: 

N a ll = N T(lS)^gg g + N x(lS)->qq + N e+e -^ qq - , (4) 

where the last piece is the most problematic background to correct for. We do this by using 
1.19 fb -1 of off-resonance data collected just below the T(AS) resonance by a process that maps 
the rj production spectra at the higher T(45) energy to the lower T(15) energy by a procedure 
analogous to what is done when one maps one normalized probability distribution into another 
one with a different domain. The remapping is done with the aid of the variable Z = E^i /Et, eam , 
the fractional energy of the rj . 

Defining three differential branching fractions dn/dZ 

dn(ggg) dB(T(15) -> ggg -> rjX) 
dZ dZ x B(T(1£) ggg) 

dn(qq) dB(T(15) -* qq -» rfX) 

dZ ~ dZ x B(T(15) -> gg) 
dn(15) _ rfg(T(15) -> VX) 
<ZZ ~ dZ 

permits us to use the Z spectrum of rj and discriminate among the various phenomenological 
models that predict different forms of H(q 2 ). Fig. 6 shows the comparison between the measured 
dn{ggg) I dZ spectrum (dots) and various phenomenological predictions (histograms). Only the 
region Z > 0.7 is relevant for comparison since the models require the rj be "fast." The model 
labeled "b" is representative of perturbative QCD 11,12 and fits the data best. There is no 
evidence for any anomalously high rj'g*g coupling at large rj energies. 
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Figure 6: Measured dn/dZ spectrum of T(1S) — ■> ggg — * rj X compared with theoretical predictions. See text for 
explanation. 
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